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q , Abstract 

►7- ' We consider a one-dimensional problem modeling two-phase flow in heterogeneous porous media 

made of two homogeneous subdomains, with discontinuous capillarity at the interface between 
them. We suppose that the capillary forces vanish inside the domains, but not on the interface. 
Under the assumption that the gravity forces and the capillary forces are oriented in opposite 
P^' directions, we show that the limit, for vanishing diffusion, is not in general the optimal entropy 

■ solution of the hyperbolic scalar conservation law as in the first paper of the series [10]. A non- 

classical shock can occur at the interface, modeling oil-trapping. 
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> ■ 1 Introduction 

The models of two-phase flows provide good first approximations to predict the motions of oil in 
the subsoil. Although the theoretical knowledge concerning the question of the existence and the 
uniqueness of the solution to such models for homogeneous porous media [4, 15] and for media with 
regular enough variations [16] is quite complete, few results are available for discontinuous media, as 
for example media made of several rock types [3, 7, 9, 11, 13, 18]. 



One says that oil-trapping occurs when some oil can not pass through interfaces between different 
rocks. Such a phenomenon plays an important role in the basin modeling, to predict the position of 
^ , eventual reservoirs where oil could be collected. As already explained in [7, 35], discontinuities of the 

!-h ' capillary pressure field can induce the so-called oil-trapping phenomenon. 

The effects of capillarity, which play a crucial role in oil-trapping, seem to play a less important 
role concerning the motion of oil in homogeneous porous media, and can sometime be neglected to 
provide the so-called Buckley-Leverett equation. 

In this paper, we show that even if the dependence of the capillary pressure with respect to the 
oil-saturation of the fluid vanishes, the capillary pressure field still plays a crucial role to determine 
the saturation profile. In order to carry out this study, we restrict our frame to the one-dimensional 
case. We will strongly use some recent results [9, 11, 13] obtained on flows in heterogeneous media 
with discontinuous capillary forces. 

We consider a one-dimensional porous medium, made of two different rocks, represented by = 
M* and Sl 2 = R+- Let tt(u,x) be the capillary pressure, then it it is well known (see e.g. the 
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introduction of the associate paper [10]) that, if both phases have different densities, the equation 
governing the two phase flow can be written 

d t u + d x (qc(u, x) + g(u, x) (1 - Cd x ir(u, x)) J = 0, (1) 

where u is the oil saturation of the fluid, q is the total flow rate, supposed to be a nonnegative constant, 
C is a constant depending on the buoyancy forces and 

c(w, x) = Ci(u), g(u,x) = gi(u), and 7r(u, x) — TTi(u) if x G Oj. 

The functions Cj are supposed to be increasing and Lipschitz continuous with q(0) = and c,(l) = 1, 
while gi are supposed to be Lipschitz continuous, strictly positive in (0, 1) satisfying g,(0) = <7,(1) = 
and 71^ are increasing Lipschitz continuous functions. 

Physical experiments suggest that the dependence of 71", with respect to u can be weak, at least 
for u far from and 1. So we want to choose 7Ti(w) = Pi, and ^(m) = Pi- The equation (1) turns 
formally to the scalar conservation law with discontinuous flux function 

d t u + d x f(u,x)=0, (2) 

where f(u,x) (resp. fi(u)) is equal to qc(u,x) — g(u,x) (resp. qci(u) — gi(u)). 

Such conservation laws have been widely studied in the last years. For a large overview on this 
topic, we refer to the introduction of [8], or in a lesser extent to the associated paper [10]. In particular, 
it has been proven by Adimurthi, Mishra and Veerappa Gowda [2] that there might exist an infinite 
number of ^-contraction semi-groups corresponding to the equation (2). Among them, in the case 
where the functions fi have at most a single extremum in (0,1), we mention the so-called optimal 
entropy solution which corresponds to the unique entropy solution in the case of a continuous flux 
function fx = fi = f. We refer to [2] and to the first part of this communication [10] for a discussion 
on the so-called optimal entropy condition. 

In the sequel of this paper, we suppose that 

(HI) for i € {1, 2}, there exists a value u* £ [0, 1) such that fi{u*) = q, fi is increasing on [0, it*] and 
fi(s) > q for all s G 1). 

We refer to Figure 1 for an illustration of the previous assumption. We denote by 

r " 

tpi(u) = C gi(s)ds. 
Jo 

For technical reasons, we have to assume that 
(H2) there exist R > 0, a > and m £ (0, 1) such that 

/io^W>? + %(l)-sr ifaebi(l)-a,¥>i(l)]. (3) 



These assumptions arc fulfilled by models widely used by the engineers, for which a classical choice of 
Ci,gi is 

_ u a * _ u a >(l-u)^ 

Cl[U, ~ u «i + |(l- u )ft' 9t[U) - 1U bu*' +a(l-u)fr' 

where o;,,/3, > 1 and a, b are given constants. 

The goal of this paper is to show that if the capillary forces at the level of the interface {x = 0} 
are oriented in the opposite sense with respect to the gravity forces (in our case Pi < P2), then a 
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Figure 1: example of functions /, satisfying Assumption (HI). Note the we have note supposed, as it 
is done in [1, 8], that fi has a single local extremum in (0, 1), but all the extrema have to be strictly 
greater than q. 



non classical stationary shock can occur at the interface. It was shown by Kaasschieter [24] that if 
the capillary pressure field is continuous at the interface (corresponding to the case P\ = P2), then 
the good notion of solution is the one of optimal entropy solution, computed by Adimurthi, Jaffre and 
Veerappa Gowda using a Godunov-type scheme [1] . We have pointed out in [10] that if the capillary 
forces and the gravity forces are oriented in the same sense, the good notion of solution is also the one 
of optimal entropy solution. If the assumptions stated above are fulfilled, if P\ < Pi and if the initial 
data Mo is large enough to ensure that both phases move in opposite directions, i.e. 

u* < u (x) < 1 a.e. in f^, (4) 

we will show that the limit is not the optimal entropy solution, but the entropy solution to the problem 

d t u + d x f t (u) = 0, 

u{x = 0~) = 1 and u(x = + ) = u%, (Piim) 
u(t = 0) = uo. 

In the sequel, we denote by a + (resp. a~) the positive (resp. negative) part of a, i.e. max(0, a) (resp. 
max(0, —a)), and for i = 1,2, for u, n G [0, 1], one denotes by 



k) = 
$. t _(u, K ) = 



fi(u)-fi( K ) if«>«, 

otherwise, 

fi(K)-fi{u) if u < K, 

otherwise, 



and 

k) = $i+(u, k) + $i-(u, k) = /,(max(«, k)) - / 4 (min(w, «)). 

We can now define the notion of solution to (Vn m ) , which is in fact an entropy in each subdomain Qi , 
with an internal boundary condition at the level of the interface. 

Definition 1.1 (solution to (Pr nn )) Let u$ £ L°°(K), u* < uq(x) < 1 a.e. in tti, A function u is 
said to be a solution of (Pn m ) if it belongs to L°°(K x R+), u* < u < 1 a.e. in Qi x (0, T), and for 



i = 1,2, for all tp G X>+(fi, x R+), for all k G [0, 1], 

(u(x, t) — n^dt'tpdxdt + / (uq(x) — n) ±, ip(x, 0)dx 

$i±(u(x,t),n)d x ip(x,i)dxdt + M fi («* - k) ± V'(0 ! t)dt > 0, (5) 
where Mf i is a Lipschitz constant of /,-, and wi = 1, U2 = 

For a given uq in L°°(R), there exists a unique solution it to (Pum) in the sense of Definition 1.1, which 
is in fact made on an apposition of two entropy solutions in R± x R + . We refer to [27, 28] and [38] for 
proofs of existence and uniqueness to solutions to the problem (Pn m ). Moreover, thanks to [12], one 
can suppose that u belongs to C(R+; L] oc {R)). 

Theorem 1.2 Let uq G L°°(M) with u* < ito < 1 a.e. m f2j, then there exists a unique solution to 
C^iim) in the sense of Definition 1.1. Furthermore, if v is another solution to (Vnm) corresponding to 
v G L°°(R) with u* < v < 1 a.e. in Hi, then for all K > 0, for all tet + 

ii pR+Mft 

{u{x, t) — v(x, t)) ± dx < / (uq(x) — wo(x)) =t da; (6) 

-ii J-R-Mft 

where Mf is a Lipschitz constant of both /j. 

Assume now that both phases move in the same direction: 

< u (x) < u* a.e. in Oj, (7) 

then it will be shown that the relevant solution u to the problem is the unique entropy solution defined 
below. 

Definition 1.3 A function u is said to be an entropy solution if it belongs to L°°(R x R+) , < u < u* 
a.e. in ftj x (0,T), and for i = 1,2, for all ip G V+(R x R + ), /or re G [0,1], 

t) — re | dtipdxdt + / |uo(#) — re| ^(x, 0)da; 



^ / ^ i {u{x 1 t),K)d x ^{x,t)dxdt + \h{K)- h{K)\j i>(0,t)dt>0, (8) 

+ie{i,2r 



n, 



Thanks to Assumption (HI), there exist no x £ [0, max it*] such that fi(x) = /2(x)> /i i s decreasing 
and / 2 is increasing on (x — <5, x + 6) for some <5 > 0. Then the notion of entropy solution described by 
(8) introduced by Towers [33, 34] is equivalent to the notion of optimal entropy solution introduced in 
[2] (see also [8]). We take advantage of this by using the very simple algebraic relation (8). 

It has been proven that the entropy solution u exists and is unique for general flux functions / [6, 
Chapters 4 and 5]. In particular, the following comparison and L 1 -contraction principle holds. 

Theorem 1.4 Let uq G L°°(R) with < uq < u* a.e. in tti, then there exists a unique entropy 
solution in the sense of Definition 1.3. Furthermore, if v is another entropy solution corresponding to 
v G L°°(R) with 0<v <u* a.e. in n u then for all R > 0, for all t £ R+ 

R rR+Mft 

(u(x,t)-v(x,t)) ± dx< / (u (x) -voix^dx (9) 

-R J-R-Mft 

where Mf is a Lipschitz constant of both f t . 
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1.1 non classical shock at the interface 



As already mentioned, the optimal entropy solution can be seen as a extension to the case of discon- 
tinuous flux functions of the usual entropy solution [25] obtained for a regular flux function. We will 
now illustrate that it is not the case with the solution to ( / Pn m )- Assume for the moment (it will be 
proved later) that in the case where Uq(x) G (u*, 1) a.e. in f2j, the corresponding solution u to (Piim) 
admits Ui as strong trace on the interface. One has the following Rankine-Hugoniot relation 

/i(wi) = f2{u 2 ) = q, 
then u is a weak solution to (2), i.e. it satisfies for all tp S T>(M. x K+): 

/ / ud t ipdxdt+ / u ip(-,0)dx + / / f(u,-)d x ipdxdt = 0. (10) 
Jm + Jr Jr Jr + Jr 



Firstly, suppose for the sake of simplicity that f\{u) = f2{u) = f{u), and that q = 0, then u* = 
for i G {1,2}. The function 

. . f 1 if x < 0, 
u(x) = < _ . r _ 
w [0 if x > 

is then a steady solution to (V\i m ) satisfying (5). However, since 

HX) ~ f( s ) < o fo r a n s g (0,1), 
1 — s 

the discontinuity at {x = 0} does not fulfill the usual Oleinik entropy condition (see e.g. [31]). This 
discontinuity is thus said to be a non-classical shock. 

Suppose now that f[(\) < and that /^(wj) > 0, then the pair (1,mJ) is a stationary under- 
compressible shock-wave, that are prohibited for optimal entropy solutions [2] as for classical entropy 
solutions in the case of regular flux functions. 

Remark. 1.5 It has been pointed out in [2] that allowing a connection (A, B), i.e. a stationary 
undercompressible wave between the left state A and the right state B at the interface lead to another 
^-contraction semi-group (see [2, 8, 21]), which is so-called entropy solution of type (A,B). However, 
we rather use the denomination non-classical shock for the connection between A and B since, as 
stressed above, the corresponding solution violates some fundamental properties of the classical entropy 
solutions. 



1.2 oil-trapping modeled by the non-classical shock 

In this section, we assume that q = 0. Let u be the solution of the problem (Pnm) corresponding to 
the initial data uq. Assume that u admits strong traces on the interface. The flow-rate of oil going 
from Qi to Q2 through the interface is given by 

/i(wi) = / 2 (u 2 ) = 0. 

Thus the oil cannot overcome the interface from fii to f^, thus if one supposes that uq belongs to 
I/°°(R), with < uq < 1 a.e., then the quantity of oil standing between x = —R (R is an arbitrary 
positive number) and x = can only grow. 

Indeed, let t-2 > t\ > 0, let ( n (x) = min(l, n(x + R) + , nx~) and 9 m {t) = min(l, m(t — m(t2 — t)). 
Choosing tp(x, t) = ( n {x)8 m (t) in (10) for m, n G N yields, using the positivity of fi 

[ ([ u(x,t)( n (x)dx) d t 8 m (t)dt+ [ 6 m (t)l- f ftiu&tydx) dt<0. 
Jt x \J-r ) Jt x \nj-i/n J 
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Since u admits a strong trace on the interface, 



1 '° 



n— *oo XI 



lim - / fi(u(x,t))dx = /i(tti) = 0. 



1/rc, 



Then we obtain 

r*2 / /-o 



t)dx d t m {t)dt < 0. (11) 
*i V-fl / 

The solution u belong to C(R+; L X (K)) thanks to [12], thus taking the limit as m — > oo in (11) provides 

o M 

u(x,ti)dx< / u(x,t2)dx. 



R J-R 



Suppose now that q > 0. Thanks to what follows, we are able to solve the Riemann problem at 
the interface for any initial data 



u (x) 



U£ if x < 0, 
u r if x > 0. 



The study of the Riemann problem is carried out in Section 5, leading to the following result. 

• If ui > u*, then U\ = 1 and u 2 = u^. We obtain the expected non-classical shock at the interface. 

• If ug < u*, then u\ = U£ and u 2 is the unique value of [0, ut,] such that 72(^2) = fi{ue)- 

Using Assumption (HI), this particularly implies that in both cases, the flux at the interface is given 
by 

fi(ui) = / a (u a ) =Gi(u e ,l) (12) 
where G\ is the Godunov solver corresponding to the flux function fx - . 



Gi(a,6) 



min fi(s) if a <b, 

se[a,b] 

max fi(s) if a > b. 

s£[6,o] 



This particularly yields that for any initial data uq G L°°(R) with < uq < 1, the restriction u\ sii of 
the solution u to ill is the unique entropy solution to 

d t u + d x f\{u) = in Qi x R + , 

u(-,0) = Uo in fix, (13) 

u(0,-)=7 inM + 

for 7=1. Since the solution u to the problem (13) is a non-decreasing function of the prescribed trace 
7 on {x = 0}, we can claim as in [10] that 

u\ n = sup { v solution to (13) } . 

7 eL°°(R + ) 

0<7<1 

In particular, u is the unique weak solution (i.e. satisfying (10)) that is entropic in each subdomain 
and that minimizes the flux through the interface. 

1.3 organization of the paper 

We will introduce a family of approximate problems in Section 2, which takes into account the capil- 
larity, with small dcpcndance e of the capillary pressure with respect to the saturation. We use the 
transmission conditions introduced in [9, 11, 13, 30] to connect the capillary pressure at the interface. 
For e > 0, the problem (V e ) admits a unique solution u e thanks to [11] and it is recalled that a 
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comparison principle holds for the solutions of the approximate problem (P e ). Particular sub- and 
super-solution are derived in order to show that if uq(x) > u J a.e. in Qi, then the limit u of the 
approximate solutions {u E ) £>0 as e tends to 0. An energy estimate is also derived. 

In Section 3, letting e tend to 0, since no strong pre-compactness can be derived on (u e ) £ > 
in Lj oc (M. x R + ) from the available estimates, we use the notion of process solution [20], which is 
equivalent to the notion of measure valued solution introduced by DiPerna [17] (see also [27, 32]). The 
uniqueness of such a process solution allows us to claim that (u £ ) converges strongly in L i 1 oc ((IR x R + ) 
towards the unique solution to (Pn m ). 

In Section 4, it is shown that if both phases move in the same direction, that is if < uq < u* 
a.e. in fli, then (u £ ) converges towards the unique entropy solution to the problem in the sense of 
Definition 1.3. 

In Section 5, we complete the study of the Riemann problem at the interface. 




2 The approximate problem 

In this section, we take into account the effects of the capillarity, supposing that they are small. We 
will so build an approximate problem (P e ), whose unknown u £ will depend on a small parameter e 
representing the dependance of the capillary pressure with respect to the saturation. We assume for 
the sake of simplicity that the capillary pressure in is given by: 

Trf = P+eu E . (14) 

It has been shown simultaneously in [9] and in [13] that a good way to connect the capillary pressures 
at the interface is to require 

^K)n7^ 2 )^0, (15) 

where u\ and u\ are the traces of u e on the interface, and where 7rf is the monotonous graph given by 

if se (o,i), 

n s {s) = { (-oo, Pi] ifs = 0, 

if s = 1. 

We suppose that the capillary force is oriented in the sense of decreasing x, i.e. Pi < P2 (the 
capillary force goes from the high capillary pressure to the low capillary pressure). Since e is assumed 
to be a small parameter, we can suppose that < e < P2 — Pi, so that the relation (15) turns to 

u\ = 1 or uf = 0. (16) 

The flux function in £li is then given by: 

F?{x,t) = Mu e )(x,t) ~ £d x Vi(u e )(x,t). 

Because of the conservation of mass, we require the continuity of the flux functions at the interface. 
Thus the approximate problem becomes 

d tU c + d x Ff = 0, 

u E {x = 0") = 1 or u e {x = 0+) = 0, , . 

F{ (0-) = Fi(0+), V ] 

U{t = 0) = Uq. 

We are not able to prove the uniqueness of a weak solution of (Vum) if the flux Ff "only" belongs 
to L 2 (fli x R+), and we will define the notion of prepared initial data, so that the flux belongs to 
L°°(Qi x M+). In this latter case, the uniqueness holds. 
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2.1 bounded flux solutions 

Wc define now the notion of bounded flux solution, that was introduced in this framework in [11, 13]. 

Definition 2.1 (bounded flux solution to (V s )) Let uq G L°°(R), < uq < 1, a function u £ is 
said to be a bounded flux solution if 

1. u e G L°°(R x R + ), < u < 1; 

2. d x ^{u £ ) G x R+) n Ll c (R + ;L 2 (iii)); 

3. u\(t) (1 — u e 2 (t)) = for almost all t > 0, where u £ denotes the trace of uf on {x = 0}. 

4. Vip e V(R x E + ), 



u £ (x,t)dtip(x 1 t)dxdt + J uo(x)ip(x,0)dx 
+ 1 V / [fi{u E )-ed x <p i {u e ))d x ^{x,t)dxdt = Q. (17) 



+ i£{l,2} 



Remark. 2.2 Such a bounded- flux u £ solution belongs to C(R+; i 1 1 oc (]R)), in the sense that there exists 
u £ in C(M + ;Lj l oc (R)) such that u £ (t) = u £ (t) for almost all t > (see [12]). More precisely, all t > 
is a Lebesgue point for u £ . So, the slight abuse of notation consisting in considering u £ {t) for all t > 
will not lead to any confusion. 

Proposition 2.3 Let u and v be two bounded-flux solutions associated to initial data uq,vq, then for 
all ^eP + (Mx R + ), 

(u — v) dtipdxdt + / (uq — Vq) ip(-, 0)dx 



I Y. I (®i±( u ' v ) - ed * (<Pi( u ) - tpiWf) d x ^dxdt > 0. 



(18) 



We state now a theorem which is a generalization in the case of unbounded domains of Theorem 3.1 
and Theorem 4.1 stated in [11]. 

Theorem 2.4 (existence— uniqueness for bounded flux solutions) Let uq g L°°(R) with < 
uq < 1 such that: 

• there exists a function u G i°°(R), with < u < 1 a.e. in R, satisfying d x ii G L 1 (R) n L°°(R) 
and such that (uq — u) G i 1 (R) 

• d x ipi(u ) G L°°(Qi) 

• lim x y u (x) = 1 or lim^^o u (x) = 0. 

Then there exists a unique bounded flux solution u £ to the problem (V £ ) in the sense of Definition 2.1 
satisfying (u £ — it) G L X (R). Furthermore, ifu £ ,v £ are two bounded flux solutions associated to initial 
data uq, vq then 

u (x) > vq(x) a.e. mR u £ (x, t) > v £ (x, t) a.e. in R for all t > 0. (19) 



Obviously, the existence of a bounded flux solution can not be extended to any initial data in L 1 (M). 
Indeed, the initial data uo has at least to involve bounded initial flux, i.e. d x ipi(uo) G L°°(M.). An 
additional natural assumption is needed to ensure the existence of such a bounded flux solution : the 
connection in the graphical sense of the capillary pressures at the interface. 

If (uq — u) and (vq — u) belong to L 1 for the same u, then (18) yields that the bounded flux solutions 
u e and v e corresponding to u$ and vq satisfy the following contraction principle: Vt £ R+, 



(u e (x,t) ~ v e (x,t)) dx < / (uq(x) — vq(x)) dx, 
i Jr 

providing the uniqueness result stated in Theorem 2.4. 
2.2 particular sub- and super-solutions 

We will study particular steady states of the approximate problem (fP e ). We will consider steady 
bounded flux solutions s e corresponding to a zero water flow rate, i.e. 

/ 4 (s £ )-e-^(s e ) =q in (20) 

For e > 0, there arc infinitely many solutions s £ of the equation (20). We will construct some particular 
solutions, that will permit us to show that the limit u as e tends to of bounded flux solutions u e 
corresponding to large initial data admits the expected strong traces on the interface {x = 0}. 

We will introduce now particular solutions of the ordinary differential equation 

y' = f l o V -\y)-q. (21) 



Lemma 2.5 Let 77 > 0, there exists a solution y n to (21) for i = 1 which is nondecreasing on (— oo, —0] 
equal to <^i(l) on [ — 77, 0) , satisfying y v (x) < <pi(l) if x < —rj and lim^^-oo y v (x) = u*. 

Proof: Consider the problem 

w'{x) = R(ipi(l) - w(x)) m if x<-rj, . , 

w(-n) = ^(1), W 

where R and m are constants given by Assumption (H2). The function 

vfi{x) = tp^l) - (R(l - m){-x - rf))^ 

is a solution of (22). Because of (H2), there exists a neighborhood (—77 — S, —77] of rj such that w ri is 
a super-solution of the problem 

y'(x) = fi°<fi 1 (y)-q if x<-r), ^ 
y(-r)) = <pi(l). 

Then there exists y v solution to (23) such that y v (x) = if x € ( — 77, 0) and 

y v {x) < w v (x) on (—1] — 5, —77]. 

In particular, y v is not constant equal to 1. Thanks to (HI), the function y v is increasing on the set 
{ x € fli I y v (x) £ (tp 1 (u*),ip }. Assume that there exists x+ < —77 such that y v (x+) = ipi{u\), 
then one sets y v (x) = •£i(u\) for all x £ (—00, If y v (x) > <pi(u\) for all x < 0, then y v is increasing 
on (—00, —77). Thus it admits a limit as x tends to —00, and it is clear that the only possible limit is 
ut □ 
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Lemma 2.6 Let 77 > 0, then there exists a solution z v to (21) for i = 2 which is nondecreasing on 
R satisfying z n {x) < <p<i ( 1 ^ 2 " 2 j if x < r), Z r, (x) > (^2 f 1 ^ 2 " 2 J if x > n and lim^^oo = ^2(1); 

lim^^.oo z(x) = u£. 



Proof: The problem 




¥>2 



admits a (unique) solution in C 1 (R, [0, 1]). Since is a constant solution of (21) for i = 2, then one 
has z(a;) > in M. Thanks to (HI), the function z is nondecreasing. This implies that it admits 
limits respectively in —00 and in +00. The only possible values for this limits are respectively itj and 
1. □ 

Proposition 2.7 Let 77 > 0, then there exists two families of steady bounded flux solutions (s £ ' I ') £>0 
and (s £,v ) e> Q tending in L\ oc {R) as s — > respectively towards 

1 if x < -77, 
if x G (-77,0), 

2 if x > 0, 

L i/ x < 0, 
1*2 if x £ (0, 77), 
L i/ x > rj. 



and 



Proof: We set 
and 

r '" (a;) = { ^(^ + 77) if x > oi (25) 

where the functions t/ 1 ' and z n have been defined in Lemmas 2.5 and 2.6. Since the functions (pi(s £,ri ) 
and c/?i(s e ' T) ) are monotone in fi{, there derivatives ^f^iGs 5 '' 7 ) an( i 7r ! Pi(s e ' v ) belong to L 1 (K), and also 
to L°°(M) because s e,T ' and s £,?? are solutions to (20). Thus they belong to L 2 (K). Hence, for fixed e, 
s £ ' n and s £,ri are bounded flux solutions to the problem (V e ). The convergence as e — > towards the 
functions s v and s n is a direct consequence of Lemmas 2.5 and 2.6. □ 

2.3 a ^((O.T);^ 1 ^)) estimate 

Our goal is now to derive an estimate which ensures that the effects of capillarity vanish almost 
everywhere in f2; x R + as e tends to 0. 

Proposition 2.8 Let u G L°°(R) with < u < 1 a.e. satisfying the assumptions of Theorem 2.4 
and let u e be the corresponding bounded flux solution. Then for all e G (0, 1), for all T > 0, there 
exists C depending only on Uo, gi,(fi,T such that 

V~e\\d x ^{u e )\\ L 2 [n ^ T)) < C. (26) 

This particularly ensures that 

£\\d x <pi(u £ )\\ L 2( aiX{0iT)) -> ase^O. (27) 
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The idea of the proof of Proposition 2.8 is formally to choose (u £ — u)ip as test function in (17) for a 
function x t— ► ip{ x ) compactly supported in fij. Using the fact that the flux Ff is uniformly bounded 
in L°°(VLi x (0,T)), we can let V tend towards XUn with Xsii( x ) = 1 if a: € Oj and otherwise, 
and the estimate (26) follows. To obtain (27), it suffices to multiply (26) by %Je. We refer to [10, 
Proposition 2.3] for a more details on the proof of Proposition 2.8. 

2.4 approximation of the initial data 

In order to ensure that the limit u of the approximate solutions u e as e — > admits the expected 
strong traces on the interface {x = 0}, we will perturb the initial data uq. 

Lemma 2.9 Let u G L°°(M) satisfying (4), then there exists (wq ,?j ) such that 

(a) . s £,ri (x) < u (x) < ^'''(x) a.e. m R, where the functions sf' 11 ands 6 ' 71 are defined in (24)-(25), 

(b) . e |j(9 x (y5i(wQ ,, ')|j ioo j ^ < C where C depends neither on e nor on rj, 

(c) . u £ n u in Lj oc (m) as (e,rj) -► (0,0). 

Proof: Let (p n ) ngN * be a sequence of mollifiers, then p„ * uq is a smooth function tending uo as 
n — > oo. Then, for e > 0, we choose n E W such that 

max jn, \\d x ipi(u * Pn)l! L -(no} - J' ( 28 ) 

and we define 

u^ix) = max{s e ' ,) (a;),min{s <r ' I '(x),Mo * p n (x)}} • (29) 
The point (a) is a direct consequence of (29). Letting (e, 77) — > (0,0) in (29) yields 

lim u e Q ' ri (x) — max {u*, min {1, uq(x)}} . 

(e,7j)->(0,0) 

Since ito is supposed to satisfy (4), this provides 

lim Un' r '(x) = un(x) a.e. in R. 

(s,ijH(0,0) 

The point (c) follows. In order to establish (b), it suffices to note that there exist an open subset ui of 
R such that u Q ,TI (x) is equal to uq * p n (x) for x £ u), and such that ^^{x) is either equal to s s ' v (x) or 
to t*>i(x) on u c = R \ cj. It follows from (28) that 

e||3»^(«o ,1, )ll£o.(n 4 n«) ^ L 

One has 

/i(uo' ,) )(x) - ed x Lp i {u e ' r, )(x) = q a.e. in S],nw c , 

thus 

£||d :E <Muo''')|| i oc (anw c ) < ||9-/i|U««,i)- 
This concludes the proof of Lemma 2.9. □ 

Definition 2.10 A function uq is said to be a prepared initial data if it satisfies (1 — uq) € L 1 (M), 
d x ipi(u ) G L°°(Sli) and 

s e ' v (x) < u (x) < s £ ' v {x) a.e. in R (30) 

for some e > 0, r\ > 0. 
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Since the function (e, rj) i— ► s £,ri is decreasing with respect to both arguments and since the function 
(£,77) 1— > s 6 ' 71 is increasing with respect to both arguments, if uo satisfies (30) for e = £q and 77 = 770, 
then uo satisfies (30) for all (e, rj) such that e < £0 and 77 < 770. So the following Proposition is a direct 
consequence from (19). 

Proposition 2.11 Let uo be a prepared initial data satisfying (30) for s — £0 and rj — r)o, then for all 
£ < £0, for all rj < rjQ, the solution u £ to (V £ ) satisfies 

sf^ix) < u £ (x,t) < s e ' v (x) for a.e. (x,t) elxl + . 



3 Convergence 

3.1 a compactness result 

Since (u £ ) e is uniformly bounded between and 1, there exists u £ L°°(R x (0,T)) such that u £ — > u 
is the L°° weak-star sense. This is of course insufficient to pass in the limit in the nonlinear terms. 
Either greater estimates are needed, like for example a _B V-cstimatc introduced in the work of Vol'pcrt 
[37] and in [10], or we have to use a weaker compactness result. This idea motivates the introduction 
of Young measures as in the papers of DiPerna [17] and Szcpessy [32], or equivalently the notion of 
nonlinear weak star convergence, introduced in [19] and [20], which leads to the notion of process 
solution given in Definition 3.2. 

Theorem 3.1 (Nonlinear weak star convergence) Let Q be a Borelian subset ofM. k , and (u n ) be 
a bounded sequence in L°°{Q). Then there exists u £ L°°(Q x (0, 1)), such that up to a subsequence, 
u n tends to u "in the non linear weak star sense" as n — > 00, i.e.: Vg £ C(R,R), 

j(«„) — > / g(u(-,a))da for the weak star topology of L°°(Q) as n — > 00. 
^0 



We refer to [17] and [20] for the proof of Theorem 3.1. 



3.2 convergence towards a process solution 

Because of the lack of compactness, we have to introduce the notion of process solution, inspired from 
the notion of measure valued solution introduced by DiPerna [17]. 

Definition 3.2 (process solution to ("Pii m )) A function u £ L°°(M. x R + x (0, 1)) is said to be a 
process solution to (Vn m ) if < u < 1 and for i = 1,2, V"0 £ T> + (f2j x R+), Vk £ [0, 1], 

(u(x,t, a) — n) ± dt')p(x 1 t)dadxdt + / (uo(x) — K^ip^x, 0)dx 

/ / $ i ±(u(x,t,a),K)d x 'ip(x,t)dadxdt + Mfc (Uj - k) ± ?/'(0, t)dt > 0, 
t + Jili Jo Jr + 

where Mf t is any Lipschitz constant of fi, U\ = 1 and U2 = u\. 



Hi Jo 



Lemma 3.3 Let uq be a rj-prepared initial data in the sense of Definition 2.10 for some rj > 0, and 
let (u £ ) e be the corresponding family of approximate solutions. Then 

u e (x,t)^l for a.e. (x,t) £ (—77, 0) x M + , (31) 

u £ {x, t) -»■ u* for a.e. (x, t) £ (0, rj) x R + . (32) 
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Proof: Firstly, since uq is a ^-prepared initial data, there exists Eq > such that 

S e °- V < U Q < S^'. 

Then it follows from Proposition 2.11 that for all e G (0, so), for a.e. (x, f)elx R+ 

/'"(a;) < u e {x,t) < s e ' n {x). (33) 
This particularly shows that for all e € (0,£o), for a.e. (x,t) G (—77,0) x K + , 

u E (a:, t) = 1, 

thus (31) holds. The assertion (32) can be obtained by using Proposition 2.7 and the dominated 
convergence theorem. □ 

Proposition 3.4 (convergence towards a process solution) Let uq be a prepared initial data in 
the sense of Definition 2.10, and let (u E )^ be the corresponding family of approximate solutions. Then, 
up to an extraction, u £ converges in the nonlinear weak-star sense towards a process solution u to the 
problem (7-WJ- 

Proof: Since u e is a weak solution of (V e ), which is a non- fully degenerate parabolic problem, i.e. 
<pj~ is continuous, it follows from the work of Carrillo [14] that u e is an entropy weak solution, i.e.: 

W> e v+(a t x r+), Vk g [0,1], 

(u s (x,t) — n) :il dtip{x,t)dxdt + / (uq(x) — n) ± ->p(x,0)dx 

R + JUi Jilt 



[$i±(u E (x,t), K ) - sd x (<pi(u £ )(x,t) - ¥>i(«;)) ± ] d x ip(x,t)dxdt > 0. 



This family of inequalities is only available for non- negative functions ip compactly supported in f2j, 
and so vanishing on the interface {x = 0}. To overpass this difficulty, we use cut-off functions Xi,s- 

Let 6 > 0, we denote by Xi.s a smooth non-negative function, with Xi,s(x) = if x ^ fij, and 
Xi,s( x ) = 1 if a; G ^i, \x\ > 5. Let tp e V + (fl x K+), then VXi,<5 G x M + ) can be used as test 

function in (34). This yields 

/ (u e - K^dtipXijdxdt + / (u - K) ± ip(-,0)xi,sdx 
+ JO, JQi 

§i±(u £ , k) - ed x (ipi(u £ ) - d x ipXi,sdxdt 



fe±(u £ , k) - e«9 x (<^(<ii e ) - ^(/s))*] ipd xXl ,sdxdt > 0. (34) 



We can now let e tend to 0. Thanks to Theorem 3.1, there exists u £ L°°(M. x R + x (0, 1)) such 
that 



lim / / (u E (x,t) — k) dtip(x,t)xi,s(x)dxdt 



0.t 

1 

± 



+ JQ.i JO 



(u(x,t,a) — k) dtip{x,i)Xi,&{x)dadxdt, (35) 



lim / / $i±(u 6 (x,t),K)d x ip(x,t)xi,s(x)dxdt 
£ - >0 Jr. Jsii 



Qi JO 



1 

$i±(u(x, t, a), n)d x ip(x, t)xi,s(x)dadxdt. (36) 
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Thanks to Proposition 2.8, one has 

£d x (ipi(u e ) — tp i (n)) ± tends to a.c. in fi, x (0,T) as e — > 0, 

then 

lim / / ed x ((fi(u e )(x,t) - ^{k^Bx (4>(x,t)xi,s(x)) dxdt = 0. (37) 

JR+ JO, 

Since Uq is supposed to be a 77-prepared initial data for some r\ > 0, we can claim thanks to Lemma 3.3 
that u E (x, t) converges almost everywhere on (— ?y, 77) x R + towards Hi if 2 G fii- Since for 8 < rj small 
enough, the support of B x xi,8 is included in the set where u E converges strongly, one has 



lim 



<f>i±(u e (x, t), K)ip(x, t)8 x Xi,s(x)dxdt 



$i±(v,i, K)ip(x,t)d x Xi,s(x)dxdt. 



We let now S tend to 0. Since Xi,s( x ) tends to 1 a.e. in Oj, (35) and (36) respectively provide 



lim lim 

(5->0e->0 



(u £ (x,t) — k) dtip{x,t)xi,s{x)dxdt 



(u(x, t, a) — K) ± dt'ip{x, t)dadxdt, 



(38) 



(39) 



lim lim 

<5->o e->0 



Si, 



$i±(M e (a;, i), K)B x ip(x, t)Xi,s{x)dxdt 



$i±(w(a;, i, a), K)B x ip{x, t)dadxdt. 



One has also 

One has 
then 



i-t-0 



lim / (uq(x) — k) ip{x, 0)xi,s(x)dx = / (uq(x) — n) tJj(x, 0)dx 



\$i±(ui,K,)\ < M h (ui - K,y 



$i±(iij, K)tp(x, t)d x Xi,s{x)dxdt 



< Mj t (ui — k) 



ip(x, t) \d x Xi,s{x)\ dxdt. 



Since |<9 x Xi,<5| tends to 8 x= q in the A4(R)-weak star sense where 

(4=0, c) , , =C(o), 

\ / M(R),C (R) 

we obtain that 

$i±(u e (x,t), K)ip(x,t)d x Xi,s(x)dxdt > 

M u (Hi - k) ± / ip{0,t)dt. 



lim inf lim 

<5->o e^O 



(40) 
(41) 



(42) 



Using (37), (39), (40), (41), (42) in (34) shows that u is a process solution in the sense of Definition 3.2. 

□ 
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3.3 uniqueness of the (process) solution 



It is clear that the notion of process solution is weaker than the one of solution given in Definition 1.1. 
We state here a theorem which claims the equivalence of the two notions, i.e. any process solution is 
a solution in the sense of Definition 1.1. Furthermore, such a solution is unique, and a i^-contraction 
principle can be proven. 

Theorem 3.5 (uniqueness of the (process) solutions) There exists a unique process solution u 
to the problem (Pn m ), and furthermore this solution does not depend on a, i.e. u is a solution to the 
problem (Vn m ) in the sense of definition 1.1. Furthermore, if uq, Vq are two initial data in L°°(M) 
satisfying (4) and let u and v be two solutions associated to those initial data, then for all t G [0, T), 

pR rR+Mft 

/ (u(x,t)-v(x,t)) ± dx< / (u {x) -v (x)) ± dx. (43) 

J-R J-R-Mft 



This theorem is a consequence of [38, Theorem 2]. Let u(x, t, a) and v(x, t, (3) be two process solutions 
corresponding to initial data uq and vq. Classical Kato inequalities can be derived in each f2j x M + by 
using the doubling variable technique: G V + (ili x R+), 



l r i 



Qi JO JO 



(u(x, t, a) — v(x, t, /?)) dtip(x, t)dad(3dxdt 



+ / (u (x) — v (x)) ip(x,0)dx 

/ $i±(u(x, t, a), v(x, t, 0))d x tp(x, tjdadfidxdt > 0. 
+ JQi Jo Jo 



The treatment of the boundary condition at the interface is an adaptation to the case of process 
solution to the work of Otto summarized in [28] and detailed in [27] leading to (see [38, Lemma 2]): 
W G V+iUi X 



1 r l 



SI; JO JO 



(u(x, t, a) — v(x, t, 0)) dttp(x, t)dad/3dxdt 



+ / (u a (x) — vo{x)) ± il>(x,0)dx 

jQi 

rl f l 

/ $ i± (u(x,t, a),v(x,t, f3))d x ij(x,t)dadpdxdt > 0. (44) 

+ JQ.i JO JO 

Choosing 

{1 if \x\<R + M f s, 
R + MfS + e - \x\ 
s — — if R + M f t < \x\ < R + M f s + e 
if \x\>R + M f s + e 

if s < t and tp £ (x, s) = if s > t as test function in (44) and letting e tend to provide the expected 
Z^-contraction principle (43). 

Finally, if u and u are two process solutions associated to the same initial data uo, we obtain a 
L 1 -contraction principle of the following form: for a.e. t G R+, 

l r i 

(u(x, t, a) — u(x, t, p))^ dadfidx < 0, 



o Jo 



thus u(x, t, a) = u(x,t,/3) a.e. in 1 x R + x (0, 1) x (0, 1). Hence u does not depend on the process 
variable a. 



15 



Theorem 3.6 Let uq be a prepared initial data in the sense of Definition 2.10, and let u e be the 
corresponding solution to the approximate problem (V e ). Then u e converges to the unique solution u 
to (Piim) associated to initial data uq in the L p ((0, T); L q (WL))- sense, for all p,q £ [l,oo). 

Proof: We have seen in Proposition 3.4 that u e converges up to an extraction towards a process solu- 
tion. The family (u £ ) £ admits so a unique adherence value, which is a solution thanks to Theorem 3.5, 
thus the whole family converges towards this unique limit u. 

Let K denotes a compact subset of R x [0, T], then one has 



(u e - ufdxdt = (u £ ) dx-2 / / u £ udx + / / vTdx. 

K JJK JJK 

Since u e converges in the nonlinear weak star sense towards u, 

lim / / (u e ) 2 dx = [[ u 2 dx 



K J JK 



Moreover, u e converges in the L°° weak star topology towards u, then 



£^0 

Thus we obtain 



lim / / u e udx = u 2 dx. 
k JJk 



lim J J (u € - ufdxdt = 0. 
One concludes using the fact the \u € — u\ < 1 for all e > 0. □ 

3.4 initial data in L°°(M) 

In this section, we extend the result of Theorem 3.6 to any initial data in L°°(M) satisfying (4) thanks 
to density argument. 

Theorem 3.7 Let uq £ L°°(R) satisfying (4), and let (wo l n) JIgN * be a sequence of prepared initial data 
tending to uq in L} oc (M). Then the sequence (u n ) n of solutions to (T-W) corresponding to the sequence 
(uo,n) of initial data converges in C(K+; Lj oc (M)) towards the unique solution to (Vnm) corresponding 
to solution the initial data uq. 

Proof: First, note that for all uo £ satisfying (4), there exists a sequence (ito,n)„ 6 pj* °^ P rc P arc d 

initial data tending to uq in Lj oc (M.) thanks to Lemma 2.9. 

Thanks to (43), one has for n, m £ W, for all t€f + 

R r R+M,t 

±j„^ / f... r„\ ... r~\\± 



(u n (x, t) - u m (x, i)) dx < (uo, n (x) — uo !m (x)) dx, 

R J-R-Mft 

then (u n ) n is a Cauchy sequence in C(M+; L ; 1 oc (R)). In particular, there exists u such that 

lim u n = u in C(M+; L} oc {R)). 

n — 'oo 

It is then easy to check that u is the unique solution to (Vn m ). □ 
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4 Entropy solution for small initial data 

In this section, we suppose that the initial data uq belongs to L 1 (R), and that 

< uq < u* a.e. in fij. (45) 

This initial data can be smoothed using following lemma whose proof is almost the same as the proof 
of Lemma 2.9. 

Lemma 4.1 There exists {uq)^ >0 C L 1 (IR) such that 

• d x ipi{ul) £ L^iSk), 

• esslim^/o uf,(x) = 1, 

• lim e ^o«o = u o in Lj oc (R). 

For all e > 0, there exists a unique bounded flux solution u £ to (V s ) corresponding to Uq thanks 
to Theorem 2.4. The following theorem claims that as e tends to 0, u £ tends to the unique entropy 
solution in the sense of Definition 1.3. 

Theorem 4.2 (convergence towards the entropy solution) Letuo £ L°°(R) satisfying (45) and 
let (uq) e be a family of approximate initial data built in Lemma 4-1- Letu £ be the bounded flux solution 
to (V s ) corresponding to Uq, then u £ converges to u in Lj oc (M. x R+) as e tends to where u is the 
unique entropy solution in the sense of Definition 1.3. 

Proof: Using the technics introduced in [10, Proposition 2.8], we can show that for all A £ [0, q] there 
exists a steady solution k| to the problem ("P e ), corresponding to a constant flux 

fi(K e x ) ~ edipi{K e x ) = A, 

and such that this solution converges uniformly on each compact subset of R* as e tends to towards 

K\(x) = min {/(k, x) = A} . 

K 

Following the idea of Audusse and Perthame [5] , we will now compare the limit u of u £ as e to with 
the steady state k\. Let A £ [0,q]. Since u e and k\ are both bounded flux solutions, it follows from 
Proposition 2.3 that for all V £V+(Rx R + ), 

{u £ - k|) ± d t ipdxdt + I (ul - k|) ± tp(-,0)dx 
+ f Y, f (^i±(u e ,Kl)-ed x (ip i (u E )-^ i {K\)) + )d x i } dxdt>Q. (46) 
Choosing A = q and ip(x, t) = (T — t) + £(x) for some arbitrary T > and some £ £ 2? + (R) yields 

/ [ (u e -K £ q )+£dxdt< [ (T-t) V / ed x {ipi{u £ )- Vi {K s q )) + d x £dxdt. (47) 

JO Jn JO i=1 2 JQi 

Since u £ is bounded between and 1, it converges in the nonlinear weak star sense, thanks to Theo- 
rem 3.1 towards a function u £ L°°(R x R + x (0, 1)), with < u < 1 a.e.. Then (47) provides 

u < K q = u* a.e. in Qj x R + x (0, 1). (48) 
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Let A G [0, q], then taking the limit for e — » in (46) yields 

|u — K\\ dtipdadxdt + / |uo — ka| 0)dx 



+ JR JO 



/ / / K\)d x ipdadxdt > 0. 

JR+ - Jf2; JO 



(49) 



Suppose that u£ > it*. Let k € [0, it*], we denote by k = f 1 (/2(«)) n [0, it*]. Then choosing A = /2(t) 
in (49), and letting e tend to gives: Vk G [0, it*], G 7_?+(R x M + ), 

m|n — K\dt4>dadxdt + / |uo — ti\ip(-, 0)dx 
isii 

|u — K\dti/jdadxdt + / |uo — K\tp(-, 0)dx 
i+ Jo.2 J Jn 2 

f / $i(it, K)d x ipdx+ ( ^ 2 (u,K,)d x iJjdx\ dadt>0. (50) 
Vjqi Jn 2 J 

It follows from the work of Jose Carrillo [14] that the following entropy inequalities hold for test 
functions compactly supported in fii: Vk G [0, 1], V?/> G 2? + (f2i x R+), 

|it E — n\dt^dxdt + / |uq — k\iJ>(-, 0)dx 

($i(u s ,k) - ed x |^i(u s ) - ^(/s))) fl^daaft > 0. (51) 

Hi 

Thus letting e tend to in (51) provides: W> G V+(Q 1 x1 + ),Vk£ [0, 1], 

I u — K\dt4>dadxdt + / \uq — K\tp(-,0)dx 

Hi JO Jf2i 

1 

$i(u, K,)d x ipdadxdt > 0. (52) 

Oi Jo 

Let S > 0, and let V G P+(R x R + ), we define 

Ipljfat) = 1p(x,t)xi,5{x), $2,5 = $ - $1,5, 

where xi,s is the cut-off function introduced in section 3.2. Then using tpij as test function in (52) 
and tp2,s in (50) leads to: 

I it — K,\d x ipdadxdt + / \u — k\$(-, 0)dx 







7 / / <&i(u, K)d x xpdadxdt 
■ JUi Jo 



+ / / / (®i(u,K)-$i(u,K))i>d x xi,5dadxdt>K(K,'ip,6), (53) 

JM + JSli JO 

where lim^o 72- (k, V'i ^) = 0- Since /1 is increasing on [0,it*] and /i([ii*, 1)) C [(7,00), cither k < it*, 
or /i(k) > /i(u*). This ensures that 

*i(«,«) = |/i(«)-/i(k)|, Vu6[0,«I], VfcG[0,«*L 

This yields 

|* 1 (u,K)-$i(«,«)| = [!/!(«) -/!(«)! -|/!(u)-/i(K)|| 

< |/i(«)-/i(«)| = |/i(«)-/a(«)|. (54) 
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Taking the inequality (54) into account in (53), and letting S — > provides: 

v« e [o, ||u||oo], W> e x>+(R x e+), 



|it — n\d x tpdadxdt + / |itQ — 0)dx 

Js. 

+ f J2 I i $i{u,K)d x i>dadxdt + \h{K)- / 2 (k)| / ^(0,-)di>0. 
Jr + 4 Jn t Jq Jo 

Using the work of Florence Bachmann [6, Theorem 4.3], we can claim that u is the unique entropy 
solution to the problem. Particularly, u does not depend on a (introduced for the nonlinear weak star 
convergence). As proven in the proof of Theorem 3.6, this implies that u £ converges in Lj oc (R x R+) 
towards u. □ 

5 Resolution of the Riemann problem 

In this section, we complete the resolution of the Riemann problem at the interface {x = 0}, whose 
result has been given in section 1.2. Consider the initial data 



u (x) 



Uf, if x < 0, 
u, r if x > 0. 



We aim to determine the traces (1*1,^2) at the interface of the solution u(x,t) corresponding to uq. 
This resolution has already been performed in the following cases. 

(a) . u\ < U£ < 1 and u* 2 < u r < 1: it has been seen that U\ = 1 and u 2 = u* 2 . 

(b) . < U£ < u\ and < u r < ui,: Since u is the unique optimal entropy solution studied in [2, 24], 

then u\ = ui and U2 is the unique value in [0, u|] such that fi(ug) = f 2(^2). 

In the cases 

(c) . u* < ut < 1 and u r = 1, 

(d) . ui = u\ and u r = 1, 

it is possible to approach the solution u by bounded flux solutions u e that are constant equal to 1 in 
fi 2 x R + . Then one obtains u\ = 112 = 1 for the case (c) and u\ = u\ and — 1 for the case (d). 

The last points we have to consider are 

(c). u\ < ui < 1 and < u r < 

(f). < ug < u* and u\ < u r < 1. 

To perform the study of the two last cases (c) and (f), we need the following lemmas that can be 
proved using similar arguments than those used in [11], particularly concerning the treatment of the 
boundary condition imposed on {x = 0}. 

Lemma 5.1 Let u r € [OjitJ)- F° r a ^ e > 0, there exists a function v £ solution to the problem 
d t v e + d x (f 2 {v £ ) - ed x <p 2 {v £ )) =0 if x > 0, t > 0, 

f 2 (v E ) - ed xV2 {v £ ) = f 2 (u* 2 ) if x = 0, t > 0, (55) 

v 6 = u r if x > 0, t = 0, 

satisfying furthermore u r < v e < u 2 and d x Lp 2 {v e ) € L°°(M + x R+). 
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Lemma 5.2 Let up £ [0, u\], u r € (1*2; 1] and let u 2 be the unique value of [OjitJ] such that f2(u 2 ) = 
fi(ui). For all e > there exists a function w £ solution to the problem 

d t W + d x (f 2 (w s ) - ed xVi (vf)) =0 if x > 0, t > 0, 
/2K) - ed x ip 2 {w e ) - / 2 (« a ) = AM tf a; = 0, <> 0, (56) 
w e — u r if x > 0, t = 0, 

satisfying furthermore u 2 < w e < u r and d x (p 2 (w 6 ) S L°°(R + x R+). 

The case (e). Assume that up > u\ and u r < u^. Let (mq)^ be a family of initial data such that 
d x fi(uQ) S L°°{£li), Uq(x) = 1 for x G (— r/,0), 1*0(2;) e [u^,l] for a.c. a; e Ox, Uq(x) = u r a - e - m ^2 
and such that 

IK -Mii(ni) < 2r7- 

Then thanks to Theorem 2.4, there exists a unique bounded flux solution u £,ri to the problem (V s ) 
corresponding to the initial data Uq. It is easy to check that the solution defined in fl 2 x R+ by the 
function v e introduced in Lemma 5.1 and coinciding in £!i x R + with the unique bounded flux solution 
corresponding to the initial data 

u o( x ) = 

In particular, as e tends to 0, it follows from arguments similar to those developed in the previous 
sections that u e,ri converges in L\ JJli x R + ) towards the unique entropy solution to the problem 
problem 

d t vP + d x f x {u r ') = ifx<0, t>0, 

if x = 0, * > 0, (57) 
if x < 0, t = 0. 
if x > 0, i > 0, 

if x = 0, t > 0, (58) 
if x > 0, i = 0. 

Note that the trace condition on the interface {x = 0} in (58) is fulfilled in a strong sense since 
u r < u(x,t) < u\ a.e. in fl 2 x R + and / 2 is increasing on [it,.,?^]. 

The solution to (57) depends continuously on the initial data in Lj oc . Hence, letting rj tend to 
in (57) provides that the limit u of u n is the unique entropy solution to the problem 

9 t u + 9 x /i(w) = ifx<0, t>0, 
u=l if x = 0, t > 0, 

u ~ ue if x < 0, t = 0. 

Note that since u{ < u? < u < 1 and mm se r Ui i] /i(s) = /i(l) = <7, the trace prescribed on the interface 
{x = 0} is fulfilled in a strong sense. This particularly yields that in the case (e), the solution to the 
Riemann problem is given by 

U\ = 1, u 2 = U 2 - 

The case (f). Following the technique used in [10] and in Section 4, there exists a unique function 
u\ solution to the problem: 

AC"!) -e-j^<pi(uf) = fi(u t ) if x<0, 
itf(O) = 1 if i=0, 
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Let u e be the function defined by 



u £ (x, t) 



u\ (x) if x < 0, t > 0, 
w e (x,t) if x > 0, t > 0, 



where w £ is the function introduced in Lemma 5.2. Then u e is a bounded flux solution to the prob- 
lem (V e ) in the sense of Definition 2.1. 



One has 
and 



— ► u t inLj- oc (ni) as e -> 0, 
ui in L i 1 oc (ri2 x R + ) as e — > 



ir 

where iu is the unique solution to 

d t w + d x f 2 (w) =0 ifx>0, £ > 0, 

w = U2 = Z^ 1 ° fl (ui) if X = 0, t > 0, 
u> = u r if x > 0, t = 0. 

Since u>(x, i) € [it2, tt r ] a.e. in x M + and since min s6 [ Uo . w ] h(s) — 72(^2) = the trace w = 1*2 

is satisfied in a strong sense on {x = 0}. This yields that the solution to the Ricmann problem in the 
case (f) is given by 

ui = ue, u 2 = /2" 1 o fi(ue). 

6 Conclusion 

The model presented here shows that for two-phase flows in heterogeneous porous media with negligible 
dependance of the capillary pressure with respect to the saturation, the good notion of solution is not 
always the entropy solution presented for example in [1, 6], and particular care as to be taken with 
respect to the orientation of the gravity forces. Indeed, some non classical shock can appear at the 
discontinuities of the capillary pressure field, leading to the phenomenon of oil trapping. We stress 
the fact that the non classical shocks appearing in our case have a different origin, and a different 
behavior of those suggested in the recent paper [36] (see also [26]). Indeed, in this latter paper, this 
lack of entropy was caused by the introduction of the dynamical capillary pressure [22, 23, 29], i.e. the 
capillary pressure is supposed to depend also on d t u. In our problem, the lack of entropy comes only 
from the discontinuity of the porous medium. 

In order to conclude this paper, we just want to stress that this model of piecewise constant capillary 
pressure curves can not lead to some interesting phenomenon. Indeed, if the capillary pressure functions 
ni are such that 7Ti((0, 1)) n 7r 2 ((0, 1)) 7^ 0, it appears in [11, Section 6] (see also [7]) that some oil can 
overpass the boundary, and that only a finite quantity of oil can be definitely trapped. Moreover, this 
quantity is determined only by the capillary pressure curves and the difference between the volume 
mass of both phases, and does not depend on uq. The model presented here, with total flow-rate q 
equal to zero, do not allow this phenomenon, and all the oil present in f2i at the initial time remains 
trapped in Oi for alH > 0. 
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